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Q ' Protoplanets accrete material from their natal protostellar disks until they are 

^ ' sufficiently massive to open a gap in the face of the disk's viscosity that arises from 

f—i ■ the magneto-rotational instability (MRI). By computing the ionization structure 
within observationally well-constrained disk models, we demonstrate that poorly 

>■ ■ ionized, low viscosity "dead zones" stretch out to 12 AU within typical disks. 

^ ■ We find that planets of terrestrial mass robustly form within the dead zones 

CN ' while massive Jovian planets form beyond. Dead zones will also halt the rapid 

^ ■ migration of planets into their central stars. Finally, we argue that the gravi- 

^ ■ tational scattering of low mass planets formed in the dead zone, to larger radii 

^ ' by a rapidly accreting Jupiter beyond, can explain the distribution of planetary 

CL|' masses in our solar system. 

o ■ 

$—1 ' 

^ ' Subject headings: planetary systems: formation - solar system: formation - plan- 

^ ■ ets: general - accretion disks - MHD - stars: pre-main-sequence 
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I 1. Introduction 

Current models for planet formation suggest either that Jovian planets are formed 
through gas accretion onto cores with ~ 10 Earth masses that are themselves assembled 
out of planetesimals (Mizuno 1980; Pollack et al. 1996), or through gravitational collapse 
when protostellar disks become cold and/or dense enough to be locally gravitationally un- 
stable (Boss 1997; Mayer et al. 2002). However, independent of how their protoplanetary 
cores were formed, planets must open gaps in their natal disk before they become isolated. 
Ultimately, it is the process of gap-opening that terminates their growth or at least severely 
restricts their final masses. A planet must be sufficiently massive to tidally open a gap in 
the face of the disk's viscosity which acts to fill in a forming gap (Lin & Papaloizou 1993; 
Ida & Lin 2004). 
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In this Letter, we calculate the gap-opening masses within protostellar disk models 
that are well constrained by the observations and analyzed in Matsumura & Pudritz (2003) 
(MP03), and show that there should be a distinct break in planetary masses akin to the 
difference between terrestrial and Jovian planetary masses. The reason for this break is that 
protostellar disks are likely to contain regions of very low viscosity - dead zones - in which 
the planetary gap-opening mass is much smaller than in regions of normal viscosity. In §2, 
we lay out the theory of gap-opening masses in turbulent disks. We updated our MP03 
disk models to include the effect of a nearby OB star, additional ionization sources, dust 
grains, as well as the turbulence driven in a dead zone by active layers (§3) and find that 
terrestrial mass planets form within the dead zones, while gas giant planets form outside of 
them (§4). Finally, in §5, we apply our results to show that dead zones can halt the rapid 
inward migration of protoplanets and that scattering of low mass planets in the dead zone 
to larger disk radii by a rapidly accreting gas giant beyond, could explain the structure of 
our own solar system. 



The minimum gap-opening mass in a fully turbulent region of a protostellar disk, Mp^turb, 
is attained when the angular momentum transport rate by a disk's turbulent viscosity be- 
comes equal to that by a planetary tidal force (Lin & Papaloizou 1993). It is readily shown 
that the gap-opening mass depends only on a disk aspect ratio h/rp (h is the local pressure 
scale height of a disk and is an orbital radius of the planet) and the turbulent viscosity 
parameter ctturb- For the case that the disk's turbulence is absent, we applied the minimum 
gap-opening mass Mp^^amp obtained by Rafikov (2002). Lin & Papaloizou (1993) also cal- 
culated a gap-opening mass in an inviscid disk by assuming that the density waves shock 
and damp immediately, but this requires a rather large planetary mass - the mass inside 
the Hill radius th = h. Rafikov (2002) showed that the density waves excited by a much 
smaller planet eventually shock and lead to the gap formation. The ratio of the gap-opening 
planetary mass for a region that is fully turbulent, compared to that for an inviscid disk, 
can be written as 



in all cases relevant to our disk model; where Q — D,Cs/ (ttGE) is the Toomre parameter that 
measures the gravitational stability of the disk, Q is the Keplerian angular velocity, and E 
is the surface mass density. For standard parameters, the critical mass ratio is of the order 
Mp^turb/Mp^damp ~ 100 which is comparable to the mass difference between terrestrial and 
Jovian planets. 



2. Gap-opening masses and disk viscosity 




(1) 
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A major source of a disk's viscosity is thought to be hydromagnetic turbulence that is 
driven by the magneto- rotational instability (MRI) (Balbus & Hawley 1998). This mecha- 
nism requires good coupling between the partially ionized gas of the disk and the magnetic 
field. Poor coupling, which occurs when this gas is not ionized sufficiently, leads to the 
formation of a so-called dead zone wherein the MRI disk viscosity is effectively zero. The 
spatial extent of a dead zone has been calculated by many authors (Gammie 1996; Sano et al. 
2000; Glassgold et al. 2000; Promang et al. 2002; Matsumura & Pudritz 2003). Physically, 
the size of a dead zone can be determined by the condition that the local growth rate of the 
MRI turbulence (~ Va/^) becomes smaller than the Ohmic diffusion rate (~ ri/h^) where 
■q is the magnetic diffusivity, and Va = B j y/A^p ~ c^turb'^s Alfven speed (p is the 

mass density, B is the magnetic flux density, and Cs is the sound speed). Detailed numerical 
calculations show that this occurs when the ratio of these two growth rates - known as the 
magnetic Reynolds number: 

Reu ^ (2) 

is less than a critical value of -ReM,crit = 10^ — 10^ (Fleming et al. 2000). We use ReM,crit — 
10^ — 10^ and ctturb = 10~^ — 1 in this paper ^. Since the magnetic diffusivity 77 is inversely 
proportional to the electron fraction, the MRI turbulence tends to be absent in poorly ionized 
regions of the disk. 

In the surface layers of a disk, MRI turbulence will always be present and can drive verti- 
cal oscillations into the dead zone below. This process leads to angular momentum transport 
whenever the dead zone is not significantly denser than the active layers: Edz/Eal < 10 
(Fleming & Stone 2003), where Sdz and Eal are the surface mass density of the dead 
zone and the active layers respectively. Therefore, the "true" dead zone should satisfy both 
conditions and we calculated it using a self-consistent model for the disk structure. 



3. Disk Model 

Most stars are thought to be formed in star clusters such as the Orion nebula cluster. A 
protostellar disk in such an environment will be irradiated by nearby luminous (10^ — lO^L©) 
OB stars in the cluster. At the outer part of the protostellar disk, the disk heating by the 
UV radiation from a central T-Tauri star (TTS) is overwhelmed by the combined radiation 
field of the nearby OB stars. As the outer disk temperature increases, the disk fiares more 
strongly and therefore planets are expected to have a higher gap-opening mass. Recently, 
Robberto et al. (2002) (RBP02) improved the self-consistent model of passive protostellar 
disks of Chiang & Goldreich (1997) by taking into account the effect of an external OB 
star. They showed that the disk aspect ratio changes significantly. We used their model to 
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calculate the vertical structure of a disk around a TTS which is exposed to OB stars ^. 

Following RBP02 and without loss of generality, we assume that a disk is oriented face- 
on with respect to an O star which is located at a typical distance of 0.1 pc from the disk and 
has the luminosity of = 6 x 10^^ erg s~^ and a Stromgren radius of it^nii = 10^^ cm. In all 
of our models, we adopt the disk surface mass density at 1 AU of So = 10^ or 10^ g cm"^ 
which is equivalent to a disk mass inside 100 AU of ~ 0.01 or O.IMq respectively (Mq is the 
mass of the Sun). These disk masses are typical for T Tauri stars (Hartmann et al. 1998; 
Kitamura et al. 2002). 

The electron fraction that controls the magnetic diffusivity of disks, is determined by 
the balance between ionization and recombination. For the ionization sources, we include 
X-rays from the central star and a nearby O star, cosmic rays, radioactive elements as well as 
the thermal collisions of alkali atoms. Among those. X-rays from the central star and cosmic 
rays are the two major ionizing sources of protostellar disks (e.g. MP03). Protostellar X-rays 
may be generated within large loops of stellar magnetic field that result from reconnection 
during magnetospheric accretion, as shown in time-dependent calculations (e.g. Hayashi et al. 
1996). We assume that magnetic loops extend out to distances (r, z) — {2R^, 2R^) from the 
stellar surface and adopt a typical T Tauri star's X-ray luminosity — 10^° erg s~^ and 
temperature kT^ = 2 keV (Feigelson, private communication). We also assume that cosmic 
rays propagate along field lines and that the strength of cosmic rays is constant across the 
disk surface. Both ionization rates as well as the uniform ionization by radioactive elements: 
^RA = 6.9 X 10~^^ s~^ arc calculated as in MP03. In this paper, we also include the thermal 
ionization of alkali metals (potassium) due to heating by the central star following Fromang 
et al. (2002). They showed that there is a magnetically active zone at the innermost radii 
(r < 0.1 AU). We find, however, that the thermal ionization effect is killed if we take 
account of the recombination on grains (see below) . Our dead zone stretches from the inner 
disk radius to a few tens of AU. We also find that the X-ray ionization by a fiducial nearby 
O star with the X-ray luminosity = 10^^ erg s~^ and temperature kT^ — 2 keV is too 
weak to affect the ionization structure of the disk. 

For recombination processes, we considered the reactions among electrons, molecular 
ions, metals and grains. At disk density n > rZcrit ~ 10^^ cm~^, grains are very effective 
at reducing the charge in disks and hence increasing the diffusivity 7]. Were it not for the 
stimulation of turbulence in the body of the disk from the envelope, our detailed calculations, 
which followed the method by Umcbayashi & Nakano (1990) and Fromang et al. (2002), show 
that very extensive dead zones - out to 16 - 29 AU - are to be expected ^. However, this 
"stimulated" turbulence limits the extent of a dead zone to smaller radii, typically 12 - 25 
AU for our two fiducial disk column densities. 
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4. Results 

We show the spatial structure of disks as well as their internal dead zones in Fig. 1 
and 2 which correspond to disk column densities of Sq = 10'^ and 10*^ g cm~^ respectively. 
Three curves show the disk height Zrem where the magnetic Reynolds number Rcm reaches 
its critical value 10^ for ctturb — 0.1, 0.01, and 0.001 from inside to outside. Also shown is 
the disk height Zgurf where the surface mass density ratio of below and above it is equal to 
10; Sbeiow/^above = 10 (a dashed line). We define the dead zone as the region where z < Zrem 
and Zrem > -Zsurf - the intersection of ^rem and Zgurf marks the outer dead zone radius As 
the magnetic field becomes stronger (the parameter aturb becomes larger), the dead zone 
becomes smaller. For a standard disk (Sq = 10^ g cm~^) with {Rcm, «turb) = (10'^, 0.01), we 
find that the dead zone stretches from the inner edge of a disk to 12 AU and encompasses 
the terrestrial planet region in our solar system (0.3 — 2 AU). For a moderately heavy disk 
(So = 10"^ g cm~^), the dead zone radius becomes 25 AU. These results agree well with 
previous works (e.g. Sano et al. 2000). It may be possible that cosmic rays do not ionize 
disks because they are swept away by disk winds. In this case, disk ionization is determined 
mainly by X-rays from the central star and an external star. The dead zone radii then 
become slightly larger (14 AU for Eq = 10^ g cm~^ and 36 AU for Eq = 10^ g cm~^). 

In Fig. 3 and 4, we show the calculated gap-opening masses of planets for disks with 
Eo = 10^ and 10^ g cm~^ respectively. Both disks are gravitationally stable (we find Qmin > 3 
in our Eq = 10^ g cm~^ model) so that planet formation by gravitational instability should 
not occur in these disk models. The upper parallel horizontal lines show the gap-opening 
masses with various viscous parameters while the bottom solid horizontal line shows the 
gap-opening masses in an inviscid disk. Vertical lines with crosses indicate the extent of a 
dead zone. A thick solid line shows the fiducial gap-opening mass throughout the disk. 

In the core accretion scenario, the minimum gap-opening mass inside the dead zone is 
dictated by wave damping (the lowermost line), while outside, the minimum gap-opening 
mass is determined by the strength of the MRI viscosity and hence the strength of the 
magnetic field (the value of aturb)- The predicted planetary mass therefore makes a distinct 
jump upwards to a value determined by the value of ccturb outside the dead zone. Our 
fiducial cases, shown in Fig. 3 and 4 with a thick solid line, corresponds to ctturb = 0.01 
outside the dead zone. This ctturb value is suggested by observations of protostellar disks 
on scales ~ 10 — 100 AU (Hartmann et al. 1998; Kitamura et al. 2002). Along this hne in 
Fig. 3, we find the gap-opening mass becomes equal to an Earth mass Me at ~ 0.7 AU and 
a Jupiter mass Mj at ~ 17 AU. Our results show that Jupiter or more massive gas giant 
planets must form outside the dead zones while terrestrial and perhaps even ice giant planets 
(see below) are likely to form within them. 
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5. An integrated picture for solar system formation 

Our results show that massive planets in our fiducial disk models form beyond 12 AU. 
This further supports the need of planet migration as an explanation of the observed exosolar 
systems. The presence of a dead zone may solve a nagging problem of migration theory - 
there is no general mechanism of halting a planet's migration into the central star. The 
standard migration picture consists of two types of migration. When a protoplanet is not 
very massive, it migrates through the disk without opening a gap (type I). As a protoplanet 
gains a sufficient mass, it opens a gap in the disk and subsequently migrates with the disk 
on a viscous timescale (type II). The type I migration is roughly two orders of magnitudes 
faster than the type II migration (Ward 1997; Terquem 2003). This means that in an inviscid 
region like a dead zone, a migrating planet will be stalled as soon as it opens a gap (Chiang 
et al. 2002). Since the presence of finite dead zones is a robust feature of the protostellar 
disks, they may act as natural barriers that prevent the rapid loss of planets into their central 
stars. 

The planetary masses predicted by our disk model (and others) increase with disk radius. 
In our solar system, this is not observed - the lower mass ice giants Uranus and Neptune 
are found at larger radii (19.2 and 30.1 AU respectively). One scenario which can explain 
this mass sequence is the photoevaporation of the disk that reduces the surface mass density 
significantly beyond ~ 15 AU within ~ 10^ years (Adams et al. 2004) but this may be too 
long for currently accepted disk lifetimes. 

We suggest an alternative - the gravitational scattering of lower mass protoplanets from 
within the dead zone to much larger radii by a gas giant located just outside of it. Numerical 
experiments (Thommes et al. 1999) have shown that a rapidly accreting Jupiter can scatter 
a more slowly growing protoplanet ary core on an interior orbit as the former's Hill sphere 
expands. The scattered low mass body will ultimately circularize its orbit by dynamical 
friction at sufficiently large disk radii to be decoupled from the scatterer. This scenario has 
the distinct advantage of building the cores of ice giants much faster because it happens 
in the inner region of the disk. This process may occur naturally in our disk model. For 
example, in Fig. 3 along the curve {Rcm, c^turb) = (10"^, 0.01), we find that a 0.74Mj planet 
is formed just outside the dead zone at r = 12 AU, while the maximum mass just inside 
the dead zone is O.OSMj - which is roughly equal to the mass of an ice giant planet. The 
inner core(s) would be scattered when the disk gas becomes sufficiently tenuous so that the 
eccentricity can be excited. 

An inevitable consequence of a dead zone is that material from the well-coupled region 
beyond its outer edge will accumulate at this interface. The increasing column density of such 
an annulus may push the dead zone outward in radius, and may even become gravitationally 
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unstable. We are currently investigating these time-dependent problems. 

In conclusion, we have calculated planetary gap-opening masses in standard mass (~ 
0.01 Mq) and moderately heavy (~ 0.1 Mq) disk models which are exposed to a nearby 
O star. With widely accepted values of (-Rbm, ctturb) = (10^,0.01), we have found that the 
dead zone stretches out to 12 AU for the standard mass disk. We have shown that there is a 
distinct and model-independent range of planetary masses within the dead zone compared to 
the well-coupled zone beyond. With {Rcm, Q^turb) = (10^, 0.01), this corresponds to terrestrial 
mass planets (2.4 x 10"^ - 8.3 x 10"^ Jovian mass planets (0.74 - 7.6 Mj). The 

robust nature of our results leads us to conclude that dead zones arc typical in protostellar 
disks and may play a central role in determining the masses of planetary families as well as 
their fates. 

We thank Eric Feigelson, David Hollcnbach, Doug Lin, Edward Thommes and James 
Wadsley for stimulating discussions and an anonymous referee for a very useful report. S. 
M. is supported by a SHARCNET Graduate Fellowship while R. E. R is supported by a 
grant from the National Science and Engineering Research Council of Canada (NSERC). 
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This preprint was prepared with the A AS lAT^jX macros v5.2. 

^Recent numerical simulations by Fleming et al. (2000) defined the magnetic Reynolds number as Rcm = 

Cgh/rj and determined the critical value in the range of 10^ — 10^. In our definition of Rcm = ct^^'^Csh/r], 
(iJcM, ttturb) = (10^,0.01) and (10^,0.01) correspond to their Rcm = 10"^ and 10^ respectively. 

^We compared three different grain models which give the emissivity difference and therefore the tem- 
perature difference of up to ~ 40 %. Two of them are RBP02's models of Type I and Type II grains whose 
radius is 0.1 /um and 0.02 /xm respectively. The other is the disk model by Chiang et al. (2001) with an 
external stellar radiation. This grain model gives the lowest temperature of all. We chose an intermediate 
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Fig. 1. — Dead zones calculated for the stan- 
dard disk model with an external star. We 
used So = 10^ g cm"^, = 10^°erg s~\ 
kT^ = 2 keV, and Rcm = 10^. The up- 
permost line shows the surface disk height; 
the lowermost line shows the pressure scale 
height while three curves show the dead zone 
boundaries for a = 0.1, 0.01, and 0.001 from 
inside to outside. The dashed line is where 

^ below above 10- 
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Fig. 2. — Dead zones calculated for a heavier 
disk of So = 10^ g cm~^ with an external 
star. For the explanation of each line, see 
Fig. 1. 



model which uses the optical depth of rj^ = and the type I grains of RBP02. 

■^To be consistent with the disk model, we assumed all grains have the same radius 0.1 /im. Our electron 

fraction for a solar nebula model is about 2 orders of magnitude larger than the one obtained by Sano et al. 
(2000) because we ignore grains of charge ±2e. With their electron fraction and our choice of Rbm = 10^, 
the dead zone would be only slightly larger. 
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Fig. 3. — Gap-opening masses for a disk 
with the surface mass density, at 1 AU, of 
So = 10^ g cm-2 (or ~ O.OIM©). The 
lowermost hne shows the gap-opening mass 
for the region with no MRI viscosity while 
the upper parallel lines show the gap-opening 
masses for a different strength of magnetic 
field: solid line is for aturb = 1, dotted line 
is for 0.1, dashed-line is for 0.01, and dot- 
dashed line is for 0.001. The dead zone for 
the magnetic Reynolds number:i?eM = 10^ is 
inside a solid line with crosses, that for 10^ is 
inside a dotted line with crosses, and that for 
10^ is inside a dashed line with crosses. The 
thick solid line is our fiducial minimum gap- 
opening mass throughout the disk for MRI 
"viscosity", aturb = 0.01. 



Fig. 4. — Gap opening masses for a disk 
with the surface mass density at 1 AU of 
So = 10^ g cm-2 (or ~ O.IMq). For 
the explanation of each line, see Fig. 3. 



